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Abstract. We present 3 methods for the subtraction of 
non-cosmic and unresolved cosmic backgrounds observed 
by the Low-Energy Concentrator Spectrometer (LEGS) 
on-board BeppoSAX. Removal of these backgrounds al- 
lows a more accurate modeling of the spectral data from 
point and small-scale extended sources. At high (>|25|°) 
galactic latitudes, subtraction using a standard back- 
ground spectrum works well. At low galactic latitudes, 
or in complex regions of the X-ray sky, two alternative 
methods are presented. The first uses counts obtained 
from two semi-annuli near the outside of the LEGS field 
of view to estimate the background at the source loca- 
tion. The second method uses ROSAT Position Sensi- 
tive Proportional Gounter (PSPC) all-sky survey data 
to estimate the LEGS background spectrum for a given 
pointing position. A comparison of the results from these 
methods provides an estimate of the systematic uncer- 
tainties. For high galactic latitude fields, all 3 methods 
give 3cr confidence uncertainties of <0.9 x 10"'^ count s~^ 
(0.1-10 keV), or <1.5 x IQ-^ count s^^ (0.1-2 keV). 
These correspond to 0.1-2.0 keV fluxes of 0.7-1.8 and 0.5- 
1.1 X 10"^'^ erg cm~^ s~^ for a power-law spectrum with a 
photon index of 2 and photoelectric absorption of 3 x 10^" 
and 3 x 10^^ atom cm~^, respectively. At low galactic lat- 
itudes, or in complex regions of the X-ray sky, the uncer- 
tainties are a factor ^-^2.5 higher. 
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1. Introduction 

The Low-Energy G oncentrator Spectrometer (LEGS) 
( Parmar et al. 1997b ) is an imaging gas scintilla- 
tion proportional counter on-board the Italian-Dutch 



BeppoSAX X-ray astronomy mission. BeppoSAX 
( Boclla et al. 1997a ) was launched into a 600 km 
3?9 inchnation orbit on 1996 April 30. The LEGS 
is one of four co-aligned narrow field instruments 
(NFIs). The other NFIs are the imaging Medium- 
Energy Goncentrat or Spectrometer (MEGS; 1.3-10 keV; 
p3oclla ct al. 1997b ), the High Pressure Gas Scintil- 
lation Proportion al Gounter (HPGSPG; 5-120 keV; 
[Manzo ct al. 1997 ) a nd the Phoswich De tection System 



(PDS; 15-300 keV 
the payload includes 
2-30 keV; Jager et al. 



Frontera et al. 1997). In addition, 
cameras (WFG; 



two wide field 

1997| ) which observe in directions 
perpendicular to the NFL The MEGS background is 



discussed in Chiappetti ct al. (1998) 
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Accurate background subtraction is important for vir- 
tually all X-ray observations and especially for studies 
of extended and diffuse objects. In this paper we dis- 
cuss background subtraction techniques appropriate to 
the LEGS. Background subtraction is complicated in the 
LEGS because (1) the wings of the mirror point spread 
function distribute a significant fraction of hard X-rays 
over the entire field of view (FOV), (2) the broad detec- 
tor response spreads a fraction of low-energy X-rays over 
the entire FOV, (3) the off-axis response of the mirrors is 
complicated and single scattered X-rays from outside the 
direct FOV can reach the focal plane ( Conti et al. 1994 ), 
and (4) at low-energies the X-ray sky has a great dea l 
of spatial structure (see e.g., Snowden et al. 199^ , 1997 ). 
This spatial structure has a complex energy dependence. 
Below 0.28 keV there is a galactic plane-to-pole variation 
of a factor '^3, while at energies between 0.5 and 2.0 keV 
the North Polar Spur and the galactic bulge dominate. 
For energies >0.5 keV and 50° < In < 290° the X-ray sky 
is relatively smooth with few discrete features visible. 

2. The LEGS 

The method of operation of the LEGS is similar to that of 
conventional gas scintillation proportional counters. An X- 
ray that passes through the ultra-thin (1.25 /im) entrance 
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window and is absorbed in the 11.4 cm diameter Xe filled 
gas cell liberates a cloud of electrons. A uniform electric 
field between the entrance window, kept at —20 kV, and 
a grounded grid causes scintillation as the electrons travel 
towards the grid. The UV light from these scintillations is 
detected by a multi-anode photo- multiplier tube (PMT). 
The LEGS is also sensitive to cosmic rays since these can 
ionize the counter gas and leave tracks through the cell. 
The 2 cm diameter gas cell entrance window is supported 
by a tungsten strongback and fine grid. The strongback 
consists of a square grid of tungsten bars with a separa- 
tion of 4'. Each strongback square is further divided into 8 
by 8 pixels by a fine grid. The overall effect of the strong- 
back and fine grid is to reduce the X-ray transmission by 
between 20-40%, depending on position within the FOV. 
The LEGS includes two ^^Fe radioactive sources which 
constantly illuminate regions of the detector outside the 
sky FOV to provide monitoring of the instrument perfor- 
mance (see Fig. ||). 

The LEGS mirror system consists of 30 nested, Au 
coated mirrors with a double cone approximation to the 
Wolter I geometry. The mirror focal length is 185 cm and 
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Fig. 1. The LEGS FOV drawn to scale. The outer circle 
indicates the 37' diameter sky FOV which contains the 
standard 16' diameter source extraction region (shaded) 
and the two background semi-annuli. The positions of the 



the geometric area 124 cm^ ( ponti et al. 199 4| |199|). The ^^^ ssp^ radioactive sources are shown as crosses 



off-axis behavior of the mirrors is complicated. Gonti et 
al. (1994) demonstrate that X-rays within a cone of solid 
angle 2° can reach the focal plane after refiection off only 
one of the mirror surfaces. As part of the BeppoSAX Sci- 
ence Verification Phase, observations were performed with 
the Grab Nebula just outside the LEGS FOV. These show 
that at an offset of 60', the Grab Nebula is visible as an ex- 
tended emission region offset in the direction of the source, 
with an intensity (5.2±J-4) x 10~^ of that on-axis. The 
spectrum of the offset emission (a power-law with a pho- 
ton index, a, of 2.2 ± 0.2) is consistent with that recorded 
from the Grab Nebula on-axis {a = 2.1). 

The LEGS is only operated during satellite night 
time. Although this reduces the observing efficiency con- 
siderably, it means that that any contribution to the 
background from scattered solar X-rays, e.g., as seen 
by the ROSAT Position Sensitive Proportional Gounter 
(PSPG) when the Sun-Eart h-satellite angle is <120° (e.g.. 



Snowden fc Freyberg 1993 ), is negligible. The BeppoSAX 



observing schedule is designed to optimize the LEGS ob- 
serving efficiency. This means that, in general, only time 
constrained or Target of Opportunity observations include 
long intervals of dark Earth pointing. The LEGS is not 
operated when the satellite passes close to the South At- 
lantic Anomaly (SAA). Due to the low BeppoSAX orbital 
inclination SAA passages result in between 5 and 12 min 
of data being lost per 96 min satellite orbit. In order to 
maximize the observing efficiency, many BeppoSAX NFI 
observations are made with a Sun-satellite-target angle 
close to 120°, which is as far from the solar direction as it 
is possible to operate. 

The LEGS energy resolution is 32% full-width half- 
maximum (FWHM) at 0.28 keV and 8.8% FWHM at 



6 keV. The FOV is circular with a diameter of 37'. The 
position resolution corresponds to 90% encircled energy 
within a radius of 2.'5 at 1.5 keV. This is a factor ^4 
worse than that of the PSPG and comparable to that of 
the ASGA Gas Imaging Spectrometer (GIS). Fig. ^ shows 
the on-axis LEGS effective area. A key scientific goal of the 
LEGS is to study sources in the energy band below the in- 
strument's G edge at 0.28 keV. This typically means that 
the absorption to a source must be ^3 x 10^^ atom cm~-^. 
The prime advantage of the LEGS, compared to previous 
high throughput X-ray detectors, is its good low-energy 
spectral resolution and its low background afforded by 
the imaging characteristics. The 0.1-2.0 keV energy reso- 
lution is a factor '--^2.4 better that that of the PSPG, while 
the effective area is between a factor ~20 and 5 lower at 
0.28 keV and 1.5 keV, respectively. The energy resolution 
is similar to that of the GIS in the overlapping energy 
range, and comparable to that of the ASGA Solid-state 
Imaging Spectrometer at energies of ~0.5 keV. 



3. Data selection and processing 

In order to study the properties of the LEGS background, 
data were extracted from 350 BeppoSAX Final Observa- 
tion Tapes, corresponding to 170 individual observations. 
These data were processed using the standard LEGS pro- 
cessing pipeline described below and cleaned event lists 
for the intervals when the instrument was observing both 
the dark Earth and the sky in the nominal operating con- 
figuration extracted. Observations with <100 counts were 
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Fig. 2. A comparison of the LEGS (solid line), ASCA SIS 
(dashed hne) and ROSAT PSPC (Dashed-dotted line) ef- 
fective areas in the energy range 0.1-5.0 keV 



excluded. This left total sky and dark Earth exposures of 
2847 ks and 498.6 ks, respectively. 

Table |l| gives the contents of the LEGS telemetry for 
each detected event. The standard LECS data processing 
package (saxledas 1.8.0) performs the following steps to 
produce cleaned, linearized event lists: 

1. Conversion of RAWX and RAWY into linear coordi- 
nates referred to as DETX and DETY using the algo- 
rithm given in Parmar et al. (1997b). 

2. Determination of the overall instrument gain (which is 
temperature dependent), using the mean energy of the 
two calibration sources as a reference. 

3. Conversion of Pulse height channels to Pulse Invariant 
(PI) channels which have a fixed relation to energy. 
This conversion takes into account the overall gain and 
a FOV position dependent gain correction. 

4. Selection of intervals with nominal instrument settings 
(such as high voltage settings), and separation into 
events that are >4° from the Earth's limb (taken to 
be good on-target intervals) , and those where the dark 
Earth occults the FOV. 

5. Removal of background events using the VETO and 
BL signals. The VETO signal is the ratio between the 
light measured by the central PMT anode and the sum 
of that measured by the surrounding anodes. It there- 
fore has a high value for on-axis X-ray sources and 
lower values for events occurring further off-axis. In 
SAXLEDAS 1.8.0 a shding cell VETO window is apphcd 
to the data. The BL signal is a measure of the dura- 
tion of the light flash produced following absorption 
in the gas cell. It depends on the penetration depth of 
an event into the cell, which is energy dependent. In 
SAXLEDAS 1.8.0 events are selected if their PI and BL 
values are within a set of 62 pre-defined boxes. 



Table 1. LECS Telemetry Contents 

Parameter Range Comment 

RAWX 0-255 Event X coordinate 

RAWY 0-255 Event Y coordinate 

TIME . . . Event Time (units of 16/is) 

PHA 0-1023 Event pulse height (energy) 

VETO 0-255 Event light distribution 

BL 0-255 Event burstlength 



The standard LECS data analysis technique is to ex- 
tract counts within a radius of 8' (corresponding to 35 
RAW pixels; see Fig. |l|) from the source centroid using 
the cleaned event lists produced by SAXLEDAS. This large 
radius ensures that 95% of the 0.28 keV X-rays are in- 
cluded in the analysis. For absorbed, or very faint sources, 
extraction radii of 4' or 6' may be appropriate. The Bep- 
poSAX pointing accuracy is such that most target sources 
are located within ~2' of the nominal position of RAWX, 
RAWY = 131, 124. 



4. Standard background 

The usual method for LECS background subtraction is 
to subtract the counts obtained in the same detector re- 
gion using a standard background exposure. This con- 
sists of the sum of the blank field exposures listed in Ta- 
ble g. The LECS data for these fields were processed using 
SAXLEDAS 1.8.0. No point sources are present in the indi- 
vidual fields with 0.1-2.0 keV fluxes >1.7 x 10"^^ erg s"\ 
or 2-10 keV fluxes of >9.7 x 10"^** erg s~^. Searches were 
also made of X-ray catalogs such as the ROSAT Bright 
Source Catalog (Voges et al. 1996) to check for the pres- 
ence of contaminating point sources in, and a couple of 
degrees outside the FOV. Note that the Gal cent-2 and -3 
fields include the complex region of sky around the Draco 
nebula (e.g., |Kerp 1994| ; [Moritz ct al. 199SD . A 22.1 ks 
BeppoSAX background exposure on 1996 August 11 at 
RA= 14i^42'"27?l, Dec= -Hl9°44'10" was not included in 
the standard background since it is in the direction of the 
North Polar Spur. Examination of the LECS spectrum 
showed that it has significantly more low-energy flux than 
the fields listed in Table 13. The remaining exposure time 
is 558.6 ks. 

The standard background 0.1-10.0 keV count rate in 
the central 8' radius source extraction region is 1.0 x 
lO^'^ s~^ keV~^ arcmin^^. The dependence of the stan- 
dard background count rate on position within the FOV 
is shown in the left panel of Fig. O. This image has been 
smoothed using a Gaussian filter with a cr of 4 RAW pix- 
els. The difference in intensity between the most and least 
intense regions is a factor 1.8. The X-ray background ap- 
pears brighter close to the center of the image due primar- 
ily to mirror vignetting. Note that the LECS mirror axis is 
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offset by 3' from the center of the FOV due to an ahgnment 
error during integration. The two bright arcs are located 
radially inwards from the ^^Fe calibration sources and re- 
sult from 5.9 keV characteristic X-rays that are absorbed 
some distance from their origin. 

Fig. |j illustrates the overall shape of the standard 
background spectrum obtained in the central 8' of the 
FOV. The Non X-ray background (NXB) obtained from 
the same region of the detector is also shown (see Sect. ^. 
The difference between the two is the contribution of the 
cosmic X-ray background (CXB), spectral fits to which 
are presented in a companion paper ( Parmar et al. 199^ . 



LEGS Standard and NXB Backgrounds 



5. Non X-ray background (NXB) 

A total of 498.6 ks of NXB data was accumulated using 
LEGS dark Earth pointings. The dependence of the 0.1- 
10 keV count rate, Ct, (s~^) on time is shown in Fig. ||, 
where each data point is averaged over the entire FOV. 
The interval without data points results from successive 
gyro failures on BeppoSAX, which led to a 3 month ob- 
serving hiatus. A gradual reduction in the counting rate 
by ~15% over an interval of ~2 years is evident. This 
may be modeled as Ct = 0.1069 - 2.2 x 10"^ AT, where 
AT is the number of days since launch. Alternatively, an 
exponential fit indicates a decay constant of 11.8 years. 
The same trend is visible in data extracted from the cen- 
tral 8' radius and background semi-annuli regions (see 
Sect. |). The intensity of the primary (5.89 keV Mn K^ 
and 6.49 keV Mn K^) radiation from the ^^Fe calibration 
sources decays with a much shorter half-life of 2.73 years. 
However, the calibration sources contain small amounts of 
other radioactive materials, the decay of which could con- 
tribute to the decrease in NXB intensity. Another factor 
may be the decreasing particle background in low-Earth 
orbit as the next solar maximum is approached, due to 
the varying height of the Earth's atmosphere. Such an ef- 
fect is evident in studies of the ROSAT High Resolution 
Imager NXB ( ^nowden 1998 ). There is little evidence in 
the LEGS data for "Long Term Enhancements" or LTEs 
as seen by the PSPC QSnowdcn ct al. 1995 ). These count 
rate increases are strongest below 0.28 keV and are ap- 
proximately uniform over an orbit. They last between 15- 
30 ROSAT orbits and have peak intensities comparable to 
the low-energy CXB. 

In the central 8' radius extraction region the 0.1- 
10 keV NXB count rate is 5.2 x 10~^ s"^ keV"^ arcmin'^. 
The dependence of the NXB on position within the FOV 
is shown in the right panel of Fig. y. The difference in 
intensity between the most and least intense regions is a 
factor 3.1. In contrast to the standard background (left 
panel of Fig. y), the NXB intensity near the center of the 
FOV does not show an enhancement. Within the central 
5' radius, any such increase is <4% of the total intensity. 

Spectra of the LECS NXB are shown in Figs. |j and 
y. The overall level is approximately constant with energy 




Energy (keV) 

Fig. 4. Rebinned spectrum of the standard background 
(filled circles) and the NXB (open circles) in the central 
8' radius of the FOV 



with 3 discrete features superposed on a smooth increase 
^4 keV. These features can be modeled as narrow Gaus- 
sian emission lines at 5.12 ±0.04 keV (at 68% confidence), 
8.51 ± 0.04 keV and 10.73 ± 0.07 keV. The first feature 
shows a clear position dependence, being stronger close 
to the calibration sources. It is almost certainly produced 
by 5.9 keV characteristic X-rays from the ^^Fe calibra- 
tion sources that penetrate deeply into the detector be- 
fore being absorbed. The detected energy is lower than 
the natural energy of the events, because X-rays absorbed 
deep within the detector produce, on average, less light 
than X-rays of the same energy absorbed close to the en- 
trance window due to the different scintillation lengths 
in the driftless gas cell (Parmar et al. 1997b). The other 
2 features do not exhibit an obvious position dependence 
within the FOV and may originate from fluorescent excita- 
tion of L-shell transitions in the tun gsten window support 
structure (see Parmar et al. 1997b ). 

To illustrate the spectral changes associated with the 
decrease in count rate shown in Fig. H, Fig. g] shows the 
spectrum of the NXB before and after day 400. This shows 
that at energies ^8 keV the long term temporal evolution 
of the LECS NXB is not strongly energy dependent and 
may be simply modeled as the change in overall normal- 
ization given above. Above 8 keV, the intensity variation 
with time is less marked. This may indicate that the fluo- 
rescent line features exhibit less intensity variability than 
the rest of the spectrum. 

The geomagnetic rigidity is a measure of the minimum 
momentum required by a cosmic particle to penetrate the 
Earth's magnetic field down to the position of the satel- 
lite. Due to the almost circular, low inclination BeppoSAX 
orbit, the variation in rigidity around the BeppoSAX or- 
bit is less than is typical for low-Earth orbiting space- 
craft such as ASCA which has an orbital inclination of 
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Fig. 3. The dependence of the LEGS standard background (left panel) and NXB (right panel) on RAWX, RAWY 
position within the FOV. Gaussian smoothing filters with a ct of 4 pixels have been applied 

Table 2. LEGS observations used to create the standard background spectrum. A target name including "sec" or 
"secondary" refers to the pri me WFG target. N is the n umber of individual pointings. Ngai is the line of sight absorption 
in units of 10^" atom cm^^ (Dickey &: Lockman 1990) 



Target 



Observation 

Start Stop 

(yr mn day) (yr mn day) 



N Position (J2000) In bn Exp. 

RA Dec (°) (°) (ks) 

(h m s) (° ' ") 



N 



gal 



Polaris Region 


1996 Jul 01 


1997 May 23 


6 


02 31 42.0 


+89 15 47 


123.3 


+26.5 


69.7 


7.0 


Gal cent-2 (sec) 


1996 Aug 23 


1996 Aug 27 


1 


16 51 19.5 


+60 11 49 


89.8 


+38.1 


82.3 


2.1 


Gal cent-3 (sec) 


1996 Aug 27 


1996 Aug 31 


1 


16 35 10.7 


+59 46 30 


89.9 


+40.2 


79.9 


2.0 


SGR-A (sec) 


1996 Oct 10 


1996 Oct 12 


1 


06 12 22.7 


-60 59 04 


270.1 


-28.2 


43.3 


4.2 


Secondary Target (SI) 


1997 Mar 18 


1997 Mar 25 


1 


17 56 46.1 


+61 11 45 


90.2 


+30.1 


94.2 


3.4 


Secondary Target (S2) 


1997 Apr 13 


1997 Apr 15 


2 


18 18 20.5 


+60 58 42 


90.2 


+27.4 


26.6 


3.8 


SDG Target 


1997 Dec 13 


1997 Dec 14 


1 


23 07 53.5 


+08 50 06 


84.4 


-46.1 


13.9 


4.7 


Secondary Target (S4) 


1998 Mar 10 


1998 Mar 22 


2 


05 52 07.9 


-61 05 35 


270.0 


-30.6 


46.7 


5.3 


Secondary Target (S5) 


1998 Aug 22 


1998 Oct 01 


4 


17 52 07.4 


+61 01 01 


89.9 


+30.6 


102.2 


3.5 



31° and experiences rigidities between 6-14 GeV c^-'^. In 
the case of BeppoSAX, the rigidity typically varies from 
10 to 16 GeV c^^ around the orbit. In order to investigate 
the dependence of the NXB spectrum on rigidity, spectra 
were accumulated for intervals when the rigidity was < 
and >13 GeV c~^. The two spectra have almost identical 
shapes, with the spectrum accumulated when the rigidity 
was <13 GeV c~^ having an overall normalization 5 ± 1% 
higher. This means that for all but the shortest observa- 
tions temporal averaging will ensure that the dependence 
of the NXB on geomagnetic rigidity can be ignored. 



6. Background semi-annuli 

The same data as in Sect. || were used to create a set 
of spectra extracted from within two semi-annuli located 
17' from the center of the FOV at positions diametri- 
cally opposite to the locations of the radioactive calibra- 
tion sources (see Fig. pi). The location and sizes of the 
semi-annuli were carefully chosen to ensure that source 
"spill-over" is only significant for bright sources, while 
the correction for mirror vignetting is kept as small as 
possible. The sum of the semi-annuli geometric areas is 
equal to that of the standard 8' radius source extraction 



LEGS background subtraction 



0.160 



_ 0.140 



w 




fll 


0.120 






CO 




CC 






0.100 


=3 




O 




O 


0.080 



0.060 



0.0 



BeppoSAX LEGS 



Dark Earth Count Rate 





»i^^?,iiWi' 



200.0 



400.0 600.0 

Days Since Launch 



800.0 



1000.0 



Fig. 5. The dependence of the FOV averaged LEGS dark Earth count rate on time. The launch date was 1996 April 
30 
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Fig. 6. The dependence of the rebinned LEGS NXB spec- 
trum, averaged over the entire FOV, on time. The filled 
circles show the spectrum obtained before day 400 and 
the open circles after day 400. The lower panel shows the 
ratio of the two spectra 



region. Gounts within the semi-annuli consist of both par- 
ticle and internal background events (i.e., the NXB), sky 
background events, and any contribution from source pho- 
tons that "spill-over" into the semi-annuli. The effects of 
"spill-over" are illustrated in Fig. where the count rate 
in the semi-annuli, Gann, is plotted against source count 




0.10 1.00 10.00 100.00 

Source Count Rate (s ') 



Fig. 7. The dependence of the 0.1-10 keV background 
semi-annuli count rate on the source count rate in the 
central 8' radius extraction region. The solid line shows 
the best-fit model discussed in the text 



rate in the central 8' radius source extraction region, Cgrc- 
Extended sources such as supernova remnants and clus- 
ters of galaxies, observations where the target source is 
offset in the FOV, and pointings near the galactic center 
are excluded. The solid line shows the fit to the expected 
linear relation: Gann = 0.0212 -t- 0.00212 Gsrc- This simple 
relation ignores the the energy dependence of "spill-over" . 
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Fig. 8. Ratio of additional counts in the background scmi- 
annuli divided by source counts as a function of energy. 
The solid line shows the fifth-order polynomial fit given in 
Table | 

Table 3. Coefficients of the fifth-order polynomial fit to 
the ratio of additional counts in the background semi- 
annuli divided by source counts as a function of energy 
in keV. These coefficients can be used to correct counts 
extracted from the background semi-annuli for source con- 
tamination, ao is the constant term 



an 


Value 


an 


Value 


ao 


-2.71 X 10"^ 


"3 


-6.13 X lO"'^ 


Hi 


+4.11 X 10"^ 


a4 


-1.01 X 10"'^ 


a2 


+9.79 X 10-* 


0-5 


+9.57 X 10"^ 



It indicates that "spill-over" adds <10% to the semi-annuli 
background for on-axis source count rates <1 s~^. 

The energy dependence of "spill-over" was investigated 
by subtracting the count rates obtained with blank fields 
from each of the semi-annuli spectra and dividing the re- 
maining counts by those of the source with the standard 
background subtracted. Fig. g| shows a polynomial fit (Ta- 
ble 0) to the ratio derived in this way as a function of 
energy. There are at least two effects which produce these 
additional counts in the background semi-annuli. At low- 
energies "spill-over" is dominated by the poor detector 
spatial resolution (oc Energy~°-^), while at high-energies 
the mirror scattering wings dominate. 

The blank field exposures listed in Table were also 
used to derive the offset correction factors which are shown 
in Fig. pi Above 4 keV, where there are few sky back- 
ground counts, the correction factors were derived by ex- 
trapolating results obtained from a Crab Nebula observa- 
tion at an offset of 11'. Multiplication of the NXB sub- 
tracted spectrum in the semi-annuli by these factors gives 
the predicted background spectrum in the central 8' radius 
source extraction region when the NXB is added. These 
correction factors are not simply the mirror vignetting 
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Fig. 9. Correction factors to be applied to a sky spec- 
trum obtained from the background semi-annuli to give 
the predicted sky background in the standard 8' radius 
source extraction region 



function given in Conti ct al. (1994) since, other factors 
such as an increase in PMT noise due to the lower signal 
levels off-axis and the position dependency of singularly 
reflected X-rays may contribute. The correction factors 
do not change significantly within the range of nominal 
target locations within the FOV of ~2' (see Sect. y). Cor- 
rection factors have also been determined for extraction 
radii of 4' and 6'. These may be obtained by scaling the 
values in Fig. |9|, which are for an extraction radius of 8', 
by the ratio of extraction region areas. The ability to ac- 
curately derive these correction factors is limited by the 
low number of counts in the standard background expo- 
sure. As more LEGS background fields become available, 
these correction factors will be updated. 

7. Background subtraction techniques 

In order to study the X-ray emission from a particular 
source it is usually necessary to subtract the X-ray and 
NXB contributions from within the extraction region. In 
the case of the LECS (and MFCS) it is important to re- 
alize that strong sources outside the FOV (and within 
2° from the instrument axis) can contribute to the back- 
ground. A search should be made for the presence of such 
sources using e.g., the HEASARC's XCOLL catalog. 

In the case of the LECS, we present 3 ways of sub- 
tracting the background. An estimate of the systematic 
uncertainty associated with background subtraction may 
be obtained by comparing results obtained using the dif- 
ferent methods. 

1.1. Standard background 

This is the simplest approach since the standard back- 
ground includes the contribution from the high (>|25|°) 
galactic latitude CXB and the NXB. Since the standard 
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background has a longer accumulation time (558.6 ks) 
than individual BeppoSAX observations, this method 
does not usually add significant noise. Since the standard 
background varies with position within the FOV (Fig. 0), 
the background spectrum should be accumulated using 
the same extraction region as the source. Note that this 
technique will underestimate the low-energy background 
for pointings in directions with bright Galactic emission 
such as the North Polar Spur (see Sect. 0). 

1.2. Background semi-annuli 

This technique has the following steps: 

1. Accumulate a spectrum of the two background semi- 
annuli during the observation of interest. These can be 
selected in RAWX, RAWY using the data extraction 
program xselect as follows: 

CIRCLE (128. 0,128. 0,80.0) 
-CIRCLEC 128 . , 128 . , 62 . 85) 
-BDX(64 . , 192 . , 128 . , 128 . 0) 
-BOX( 192 . , 64 . , 128 . , 128 . 0) 

2. For source count rates ^1 s~^, subtract the NXB back- 
ground in the standard extraction region from the to- 
tal counts in the same region. Then multiply the sub- 
tracted counts by the function shown in Fig. g to derive 
the "spill-over" contribution. This is then subtracted 
from the extracted semi-annuli spectrum. 

3. Subtract the standard NXB spectrum of the semi- 
annuli. 

4. Multiply the remaining spectrum by the correction fac- 
tors shown in Fig. ||. This provides an estimate of the 
CXB in the source region. 

5. Add the standard NXB spectrum in the source region 
to this spectrum. This is then the final background to 
be subtracted before scientific analysis of the source 
properties. 

7.3. Scaled ROSAT PS PC all-sky survey background 

ROSAT PSPC all-sky survey diffuse background maps 
are available in 7 energy bands with a pixel size of 
12' X 12' covering ~98% of the sky The effects of dis- 
crete source counts, non X-ray contamination and X-rays 
of solar system origin have been eliminated to the great- 
est possible extent. The contributions of all the sources 
listed in Table have been excluded. These maps can 
be used to predict the variation with pointing position 
of the LEGS background. The PSPC count rates closest 
to the centers of each of the fields included in the LEGS 
standard background were first determined. The mean, 
LEGS exposure weighted, PSPC 0.1-2.5 keV count rate is 
1.2 X lO^'^ s^^ arcmin^^. 

In order to understand how best to model the PSPC 
data from different sky positions, 7 channel PSPC spectra 



were determined for the fields given in Tables || and |[ 
These were then fit with a number of different mod- 
els including power-laws, thermal bremsstrahlungs, bro- 
ken power-laws with the hard component a fixed at 1.4 
(the value expected above 1 keV), and power-laws with 
a = 1.4, together with thermal components. In all cases 
low-energy absorption was included. It was found that (1) 
the power-law and (2) the fixed power-law together with 
a thermal component give consistently the best fits for al- 
most all the fields. These two models are therefore used 
in the background analysis except for the Gal cent-3 field, 
where a broken power-law, with hard index a fixed at 1.4 
is preferred. 

In order to use this technique to predict the LEGS 
background spectrum at a particular pointing position, 
the following steps are necessary: 

1. Determine the ROSAT PSPC count rates in the 7 en- 
ergy bands at the required position (see Sect. ^). 

2. Fit spectral models to the derived spectrum using the 
appropriate ROSAT all-sky survey response matrix. 

3. Select the model with the lowest reduced x^- 

4. Use the standard LEGS response matrix and the pa- 
rameters of the selected model fit to calculate a pre- 
dicted background spectrum for the position of inter- 
est. 

5. Use the standard LEGS response matrix and the GXB 
fit parameters given below to calculate a predicted 
spectrum for the standard background field. 

6. Divide the predicted LEGS spectrum by that predicted 
for the standard background to derive a set of PI chan- 
nel dependent correction factors. 

7. Subtract the standard NXB spectrum (see Sect. ^ 
measured in the source extraction region from the ob- 
served standard background spectrum. 

8. Multiply the subtracted spectrum by the correction 
factors calculated above. 

9. Add the standard NXB measured in the source extrac- 
tion region from the standard background spectrum to 
the scaled spectrum. This is then the final background 
to be subtracted before scientific analysis of the source 
properties. 

The LEGS GXB spectrum is discussed in a compan- 



ion paper (Parmar et al. 1999). In the 0.1-7.0 keV energy 
range it is well modeled by an absorbed power-law to- 
gether with two solar abundance thermal bremsstrahlung 
components (the MEKAL model in xspec). When this 
model is fit to the exposure weighted PSPC spectrum 
corresponding to the standard background, the best- 
fit results are a power-law a of 1.48 with normaliza- 
tion of 1.07 X lO^'^ photon cm^^ s^^ keV^^, a hard 
bremsstrahlung component temperature (kT) of 0.80 keV 
with normalization of 1.5 x 10^^ photon cm~^, and a soft 
bremsstrahlung component kT of 0.159 keV with normal- 
ization of 1.9 X 10~^ photon cm~^. For both absorbed 
models, Nh equivalent to 3.7 x 10^° atom cm~^ is required. 
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These fit parameters should be used to generate the pre- 
dicted standard LEGS background spectrum. 



8. Background subtraction comparison 

The accuracy of the different background subtraction 
techniques described in Sect, ^was investigated using both 
blank field and exposures containing sources. The blank 
field exposures are listed in Table and the source expo- 
sures in Table 0. 



8.1. Target selection 

GAL 83 and GAL 87 are two "supersoft" sources located 
in the Large Magellanic Gloud. GAL 83 is only detected in 



the LEGS between 0.1-0.6 keV (Parmar et al. 1998) and 



GAL 87 between 0.2-1.0 keV ( [Parmar et al. 1997aD. H Z 43 
is a nearby hot white dwarf (e.g., Barstow et al. 199^ ) that 
is often used as an extreme ultra violet or soft X-ray "stan- 
dard candle" in X-ray astronomy. Due to the softness of its 
spectrum, it is only detected below 0.28 keV in the LEGS. 
These 3 targets were chosen since they are only detected 
in a narrow energy range of the LEGS, allowing the rest 
of the spectrum to be used to estimate the background 
subtraction quality. 

4U 1630—47 is a recurrent X-ray transient located close 
to the galactic plane. The BeppoSAX observation was de- 
signed to detect quiescent emission between outbursts, but 
none was found ( Oosterbroek et al. 1998 ). X 1755— 338 
is a bright X-ray dipping source which was observed 
t o have turned-off p rior to the BeppoSAX observation 
( Roberts et al. 1996|) . These two fields were chosen as tests 
of the background subtraction techniques in complex re- 
gions of the sky at low galactic latitudes. 

NGG 7172 is a Seyfert 2 galaxy, whose nuclear con- 
tinuum is seen through an Nh ~ 10^'^ atom cm^^ neu- 
tral absorber. A soft excess observed by ASGA was as- 
sociated with diffuse emissi on from the group HG G 90, in 
which NGG 7172 is located ( [Guainazzi et al. 199"8| ). In ad- 
dition, two relatively bright (~0.25 count s~-^) sources are 
included as a control sample, to show the relative inde- 
pendence of the results on the choice of the background 
subtraction method at higher flux levels. These are the 
Seyfert 2 Galaxy NGG 1068 ( puainazzi et al. 1999 1) and 
the coronal X-ray source VY Ari (Favata et al. 1997). 



8.2. Results 

In Fig. no, observed LEGS count rates in five energy bands 
(0.1-0.3, 0.3-1, 1-2, 2-5, and 5-10 keV) are shown for 
the sample sources. No correction for the instrumental re- 
sponse has been applied. All results in this section are dis- 
played in this way to allow ease of comparison since count 
rate is a quantity which can be easily estimated when 
comparing the behavior of a weak X-ray source with the 
instrumental performances and limitations. 




1.5x10 



Energy (keV) 

Fig. 12. Exposure time weighted sum of the 8 blank field 
exposure residuals when the 3 background subtraction 
techniques are applied. Empty squares indicate the stan- 
dard background, filled circles the semi-annuli, and crosses 
the ROSAT-scaled techniques 



Fig. 11 shows the residuals for 8 blank fields when each 



of the 3 background subtraction techniques are applied. 
(The shortest exposure blank field, the SDG Target, in 
Table |2| was excluded.) In the case of a good background 
subtraction, the residuals should exhibit a Poissonian dis- 
tribution centered on zero, and the standard deviation of 
the residuals provides an estimate of the systematic uncer- 
tainties associated with each technique. These results are 
summarized in Fig. |lj, where the exposure time weighted 
sum of the background field residuals, shown individually 
in Fig. |l^, are plotted. The average mean count rates, n, 
in the 0.1-10 keV energy range are (— 1.5±0.9), (3.6±1.6) 
and (— 1.3± 1.0) x 10~^ s~^, for the standard, semi-annuli 
and ROSAT-scaled technique, respectively. Tabic g lists 
the values of | /i | -I- 3cr (cr is the standard deviation) for the 
3 techniques in 3 energy ranges. These values may be inter- 
preted as 3a estimates of the systematic uncertainties of 
each the background subtraction techniques at high galac- 
tic latitudes. In the 0.1-2 keV energy range, the values in 
Table || correspond to fiuxes in the range 0.7-1.8 and 0.5- 
1.1 X 10""'^^ erg cm~^ s~^ for a power-law spectrum with 
a = 2.0 and photoelectric absorption, Nh, of 3 x lO^*' and 
3 x 10^^ atom cm^^, respectively. In the 2-10 keV band, 
for an Nh of 3 x 10^*^ atom cm~^, the limiting fluxes are 2- 
3 and 1.3-3 x 10~^^ erg cm~^ s^^ for sources with a = 2.0 
and 1.5, respectively. 

In the case of the blank fields, the standard back- 
ground provides the most consistent (i.e., lowest values 
of |/i| + 3a) subtraction. This is unsurprising given that 
the standard background itself includes the comparison 
fields (see Table |2|). In addition, the blank field expo- 
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Table 4. Observations used for background accuracy checks. CR is LEGS count rate including the contribution of the 
background in the 8' radius extraction regions 



Source 


CR 


Date 


In 


bii 


Remarks 


Name 


(S-) 


(yr mn day) 


n 


n 




CAL87 


0.031 


1996 Oct 27 


281.8 


-30.7 


Supersoft source 


4U 1630-47 Field 


0.040 


1997 Mar 26 


336.9 


+0.3 


X-ray transient in quiescence 


CAL83 


0.043 


1997 Mar 07 


278.6 


-31.3 


Supersoft source 


X 1755-338 Field 


0.045 


1998 Apr 10 


357.2 


-4.9 


X-ray binary in quiescence 


NGC7172 


0.059 


1996 Oct 14 


321.2 


-18.6 


Active Galactic Nuclei 


NGC 1068 


0.142 


1996 Dec 30 


172.1 


-51.9 


Active Galactic Nuclei 


HZ 43 


0.215 


1998 Jan 27 


54.1 


+84.2 


Hot white dwarf 


VYAri 


0.303 


1996 Sep 06 


150.6 


-25.4 


Coronal X-ray emission 



0.1 
0.01 



K 



" 10 



c^ 



-3 



-4 



10 
0.1 

^ 0.01 

u 10 ^ 

10-' 

0.1 



0.01 



" 10 



-3 



-4 



10 
0.1 

^3 0.01 

o 10"^ 

10-" 



Cal87 



HZ43 




NGC1068 




VYAri 

I I I I 



+^ 



X1755-338Field 



NGC7172 



0.1 1 

Energy (keV) 



10 



0.1 1 

Energy (keV) 



10 



Fig. 10. Count rates (CR) for the source sample in the energy ranges 0.1-0.3, 0.3-1, 1-2, 2-5, and 5-10 keV. The counts 
include contributions from the source (if present) and any extended emission components within the 8' extraction radii. 
The ordinate uncertainties are smaller than the symbol size 



sures are at high (>|25|°) galactic latitude and avoid fea- 
tures such as the North Polar Spur, and so are expected 
to be broadly similar. The residuals obtained with the 
semi-annuli method are the largest. An investigation re- 
veals that this may be in part due to uncertainties in the 



NXB spectra of the two semi-annuli which contain <8800 
counts. It is expected that this will improve as subsequent 
dark Earth exposures are added to the already existing 
data set. The ROSAT-scaled background works consis- 
tently well within the limits specified in Table y. 
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Fig. 11. Count rate residuals (CR) as a function of energy for 8 blank fields from Table g, when the 3 background 
subtracted techniques are applied. Empty squares indicate the standard background, filled circles the semi-annuli, and 
crosses the ROSAT-scaled techniques 



Fig. |lj illustrates the application of the 3 background 
subtracted techniques to the sample sources. For each 
source the 3 techniques were applied and background sub- 
tracted count rates in the 5 energy bands calculated. The 
difference with respect to the mean of the 3 techniques is 
plotted in Fig. ^. In 6 out of 8 cases no significant de- 
viations at levels higher than a few xlO^'^ count s~^ are 
present. This is comparable to the results on the blank 
sky fields and indicates that all three methods work well 
here. The differences between the three methods may be 
used to estimate the systematic uncertainties associated 
with background subtraction. 

In the case of the two fields close to the galactic 
plane (4U 1630-47 and X 1755-338), differences at a 
level ^lO^'^ count s^^ are evident. This is a factor '--^2.5 
larger than with the blank field exposures and is proba- 
bly due to incorrect estimation of the contributions of the 
hard diffuse emission associated with the Galactic ridge, 
or unresolved point-sources. In addition, in the case of 
4U 1630-47 the 50 mCrab (or ~10 LEGS count s'^) X- 



ray binary 4U 1624—49 is located 2° away. X-rays from 
4U 1624—49 that undergo a single mirror reflection may 
provide a small contribution to the 4U 1630—47 field. 

In both these cases, the standard background system- 
atically underestimates the background (as expected) , and 
the scaled ROSAT method overestimates the background 
somewhat, when compared to the mean values. The reason 
for this is unclear, but it may be partly due to incorrect 
ROSAT source subtraction in the complex fields. In the 
case of the X 1755—338 field, the residuals obtained using 
the semi-annuli method are smoothly distributed between 
those obtained using the other two methods. This is not 
the case in the 4U 1630—47 field where the residuals de- 
viate from the mean values above 2 keV. These results 
suggest that the semi-annuli gives the most reliable re- 
sults for complex fields, with the scaled ROSAT method 
being the second most reliable. These two cases illustrate 
the difficulty in obtaining a good background subtraction 
for sources located at low galactic latitudes, or in complex 
regions of the X-ray sky. 
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Fig. 13. Count rate difference compared to the mean (CR) as a function of energy when the 3 background subtraction 
techniques are applied to each of the sample sources. Empty squares indicate the standard background technique, 
filled circles the semi-annuli, and crosses the ROSAT-scaled technique. Note that the ordinate extrema are a factor 2 
larger than in Fig. [ll| 



Table 5. The values of |/i| + 3a for the 3 background sub- 
traction techniques in 3 energy ranges for the blank fields 



Method 



0.1-10 keV 0.1-2 keV 2-10 keV 



Standard 


0.8 


Semi-annuli 


1.4 


ROSAT-scaled 


0.9 



0.7 
1.9 
1.3 



1.3 
1.2 
1.3 



Finally, the effects of the long term decrease in NXB 
intensity evident in Fig. a were evaluated using all three 
proposed background subtraction techniques. In all cases 
the differences in background residuals were substantially 
smaller than the values given in Table g. This means that 
currently no time dependent corrections need to be ap- 
plied to the LEGS background. 



9. Data products 

The standard background and NXB cleaned event lists 
will be made available by the BeppoSAX Science Data 
Center in Rome. These files will allow users to generate 
appropriate background and NXB spectra. 

In addition, the Max-Planck-lnstitut fiir Extrater- 
restrische Physik, in Garching will make the 7 energy 
channel source subtracted PSPC maps available through 
the World Wide Web so that users can use the scaled 
ROSAT background subtraction technique. 
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